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Abstract. We prove that for any fixed d the generating function of the projection 
of the set of integer points in a rationai d-dimensionai poiytope can be computed in 
poiynomiai time. As a coroiiary, we deduce that various interesting sets of iattice 
points, notabiy integer semigroups and (minimai) Hiibert bases of rationai cones, 
have short rationai generating functions provided certain parameters (the dimension 
and the number of generators) are fixed. It follows then that many computational 
problems for such sets (for example, finding the number of positive integers not 
representable as a non-negative integer combination of given coprime positive integers 
ai,. . . ,aii) admit polynomial time algorithms. We also discuss a related problem of 
computing the Hilbert series of a ring generated by monomials. 

1. Introduction and Main Results 

Our main motivation is the following question which goes back to Frobenius and 
Sylvester. 

(1.1) The Frobenius Problem. Let ai,... , be positive coprime integers and 


let 



be the set of all non-negative integer combinations of ai,... , a^, or, in other words, 
the semigroup S C Z_|_ of non-negative integers generated by ai,... , a^. What does 
S look like? In particular, what is the largest integer not in S'? (It is well known 
and easy to see that all sufficiently large integers are in S). How many positive 
integers are not in S? How many positive integers within a particular interval or a 
particular arithmetic progression are not in S? 
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One of the results of our paper is that for any hxed d “many” of these and 
similar questions have “easy” solutions. For some of these questions, notably, how 
to hud the largest integer not in S, an efficient solution is already known [K92]. 
For others, for example, how to hud the number of positive integers not in S, an 
efficient solution was not previously known. 

With a subset S C we associate the generating function 

/(S;x) = 5; I™, 

mES 

Clearly, the series converges for all x such that |a:| < 1. We are interested in hnding 
a “simple” formula for f{S-,x). 

(1.2) Examples: d = 2 and d = 3. Suppose that d = 2, that is, S is generated 
by two coprime positive integers ai and 02 - It is not hard to show that 

“ (1 -a;“i)(l 

Suppose that d = 3, that is, S is generated by three coprime positive integers 
ai, 02 and as. Then there exist (not necessarily distinct) non-negative integers 
Pi,P 2 ,P 3 ,P 4 and ps, which can be computed efficiently from ai,a 2 and as, such 
that 

1 '-rPl nnP2 nnPZ _(_ _|_ rpP^ 

p , ^ .. -L kKj kKj kKj I vC/ I vC/ 

^ ( 1 -a:“i)(l-a:“ 2 )(l-a:“3) 

This interesting fact is, apparently, due to G. Denham [D96]. For example, if 
a = 23, 6 = 29 and c = 44, then (thanks to a MAPLE program written by J. 
Stembridge), pi = 161, p 2 = 203, ps = 220, p 4 = 249 and ps = 335. 

The idea of Denham’s proof is to interpret /(S'; x) as the Hilbert series of a graded 
ring M = C[t“b This ring M can be considered as a graded module over 

the polynomial ring R = C[a:i, 2 : 2 , 0 : 3 ] graded by degXi = Oi and acting on M 
by Xit"'^ = Since the projective dimension of M is 2, the Hilbert-Burch 

Theorem allows us to construct explicitly a projective resolution of M and then to 
compute the Hilbert series from it, cf. Section 20.4 of [E95]. 

We also note that a slightly weaker form of this result is obtained by elementary 
methods in [SW 86 ]. 

What happens for d = 4 (or larger)? Clearly, since S contains all sufficiently 
large numbers, f{S-,x) is a rational function of the type 

(1.3) f{S-,x) = pn{x) + - - , 

1 — a: 

where N is the largest integer not in S and Pn{x) is a polynomial of degree N. 
Can we hnd a shorter formula for /(S'; 2 ;)? 

We need some standard dehnitions from computational complexity theory (see, 
for example, [P94]). 
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(1.4) Definitions. We define the input size of an integer a as the nnmber of bits 
needed to write a, that is, ronghly, l + log 2 |a|. Hence the inpnt size of the seqnence 
ai,... ,ad will be ronghly d + log 2 aj. We are interested in the complexity of 
an algorithm which compntes /(S'; x) from the inpnt ai,... , a^. The algorithm is 
called polynomial time provided its rnnning time is bonnded by a certain polynomial 
in the inpnt size. 


We show that for any fixed d there is a mnch shorter formnla for f{S',x) than 
that given by (1.3). 


(1.5) Theorem. Let us fix d. Then there exists a positive integer s = s{d) and a 
polynomial time algorithm, which, given the input ai,... ,ad, computes f{S-,x) in 
the form 


f{S;x) = 


x^ 


Oti 


iei 


(1 — x^ii) ■ ■ ■ (1 — X^is) ’ 


where I is a set of indices, ai are rational numbers, pi and btj are integers and 
bij 7 ^ 0 for all i, j 

In particnlar, the nnmber |/| of fractions is bonnded by a certain polynomial po/y 
in the inpnt size, that is, in d + The degree of poly and the nnmber 

s = s{d) both grow fast with d, ronghly as d^^^\ However, for any fixed d, the 
formnla of Theorem 1.5 is mnch shorter than that of (1.3), in fact, exponentially 
shorter. Indeed, by [EG72] it follows that for any fixed d, the integer N in (1.3) 
can be as large as 0{t^), where t = maxjai,... , ad}. Thns the length of formnla 
(1.3) is qnadratic in t, that is, exponential in the inpnt size. For d = 4, there 
are examples (see [SW86]) showing that if the denominator of f{S',x) is chosen 
in the form (1 — x°'^){l — x‘^^){l — a:“®)(l — x‘^'^) then the nnmber of monomials 
in the nnmerator can grow as fast as \/i for t = minjai, a 2 , as, a 4 }, which is still 
exponential in the inpnt size. 

Theorem 1.5 is a special case of a more general resnlt. Let S' C be a (finite) set 
of integer points. For an integer vector m = (pi,... , pd) and (complex) variables 
X = (xi ,... , Xd), X G C'^, let 




■ ■ ■ X 


Ud 

d 


denote the corresponding monomial. We let x^ — 1. Let ns consider the Lament 
polynomial 

/(S;x)= 

mES 


This a priori “long” polynomial can sometimes be written as a “short” rational 
fnnction 


= Y 


fPi 


ai 


iei 


(1 


-)...(! 
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where G Q, Pi, bij G and hij 7 ^ 0 for all i,j. The motivating example is the 
set S = { 0 , 1 , 2 ,... , n}, for which we have 


fc =0 


1 - 
1 — X 


Thns, for this particnlar S', the long polynomial f{S-,x) can be written as a short 
rational fnnction in x. Indeed, writing /(S; x) as a polynomial reqnires, ronghly, 
0(n logn) bits whereas writing f{S',x) as a rational fnnction reqnires only O(logn) 
bits. A more general example is given by the set of integer points in a rational 
polyhedron. 

(1.6) Definition. Let Cl, ... , Cn G be integer vectors and let /li,... , /9n G Z 
be integers. The set 

P = |a: G : (c^, x) < (3i for i = 1,... , n| 

is called the rational polyhedron dehned by {ci, (3i}. Again, we dehne the inpnt size 
of P as the nnmber of bits needed to dehne P. That is, if Ci = ( 7 ^ 1 ,... , 7 ^^) then 
the inpnt size of P is ronghly 


n n d 

nd + Z log2 \/3i\ + ZZ log2 hij\- 

i=l i=l j=l 

A bonnded rational polyhedron is called a rational polytope. 

In [BP99] it is proved that for any hxed d, if P C is a rational polyhedron 
which contains no straight lines then for A = P fl Z^ the expression 

/(S;x)= 5^ x™ 

mePnzd- 

can be written as a short rational fnnction. We give the precise statement in 
Theorem 3.1. 

The main resnlt of this paper is that the projection of the set of integer points 
in a rational polytope has a short generating fnnction as well. More precisely, let 
T : be a linear transformation snch that T(Z'^) C Z^. Thns the matrix 

of T (which we also denote by T) with respect to the standard bases of R'^ and R^ 
is integral. The inpnt size of T is dehned similarly as the nnmber of bits needed to 
write T. Thns, if T = (tij): i = I,... , /c and j = I,... , d then the inpnt size of T 
is ronghly kd + J2i=i log 2 l^bl- Let S = T{P fl Z'^), S C Z^, be the image of 
the set of integer points in P. We prove the following resnlt. 
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(1.7) Theorem. Let us fix d. There exists a number s = s{d) and a polynomial 
time algorithm, which, given a rational polytope P and a linear transformation 
T : Mfi such that T(Z^) C , computes the function /(5';x) for S = 

T{P n Z‘^), S <zZ^ in the form 




fPi 


CXi 


iei 


(1 — x“ii) ■ ■ ■ (1 — X“is) ’ 


where ai G Q, pi, aij G and aij 7^ 0 for all i,j. 

In particular, the number |/| of fractions in the representation of /(5';x) is 
bounded by a certain polynomial in the input size of P and T. We do not discuss 
the exact dependence of s{d) on d but note that a rough estimate suggests that s 
can be chosen about d^^^\ 

We obtain Theorem 1.5 as a simple corollary of Theorem 1.7 (see Section 6). In 
Section 7, we discuss other interesting sets which possess short rational generating 
functions, such as the (minimal) Hilbert bases of rational cones and “test sets” in 
parametric integer programming. We also discuss a related problem of finding a 
short formula for the Hilbert series of a ring generated by monomials. 

What can we do with rational generating functions? As is discussed in Section 3, 
we can efficiently perform Boolean operations on sets given by their short rational 
generating functions. In particular, if 81,82 C Z*^ are two sets of integer points 
given by their generating functions /(S'i;x) and /(S' 2 ;x), we can compute the 
generating functions f{8i fl 5'2;x), f{8i U 5'2;x) and f{8i \ 5'2;x) in polynomial 
time (see Theorem 3.6). Also, by specializing at x = (1,... ,1), we can count 
points in polynomial time in finite sets given by their generating functions (this is 
not immediate since x = (!,...,!) is a pole of each fraction in the representation 
of /(S';x), cf. Theorem 2.6). 

Let f{8',x) be the generating function of Theorem 1.5. Then, for the complement 
8 = Z_|_ \ 8, we compute the generating function f{8-,x) = (1 — a:)“^ — /(S'; x) and 
then compute the number of non-negative integers not in 8 by specializing f{8',x) 
aX X = 1. Given an interval [a, b] C Z_|_, for S' = Sfl [a, b], we can compute /(S'; x), 
and, specializing a: = 1, we can obtain the number of points in S inside the interval 
[a,b]. 

The proof of Theorem 1.7 combines several methods. First, it uses some tech¬ 
niques of working with short rational generating functions developed by the first 
author, see [BP99] and Sections 2 and 3. Second, it uses some “flatness”-type ar¬ 
guments from the geometry of numbers, see, for example, [GLS93] and Section 4. 
Finally, it relies on parametric integer programming arguments developed by R. 
Kannan, L. Lovasz and H. Scarf, see [K92], [KLS90] and Section 5. The crucial 
step of bringing the three ideas together and obtaining the proof of Theorem 1.7 is 
made by the second author (Section 6). 


Remark. When a lemma or a theorem states that “there exists a polynomial time 
algorithm”, the actual algorithm is either provided in the proof or a suitable refer¬ 
ence is given. 
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2. Rational Functions and Monomial Substitutions 

In this section, we develop certain methods of specializing rational fnnctions 
/(x), X G of the type 


/(x) = 


E 

iei 


rVi 


aj 


(1 — x ““) ■ ■ ■ (1 — ) ’ 


where / is a hnite set of indices, ctj G Q, Pi, aij G and ^ 0 for all i, j. We 
hx an npper bonnd k > k{i) on the nnmber of binomials in every denominator bnt 
allow the nnmber of variables d, the nnmber |/| of terms, the coefficients cti and 
the vectors Pi,aij to vary. Moreover, to simplify the notation somewhat, we will 
consider the case of all k{i) being eqnal to a nnmber /c, so 


( 2 . 1 ) 


/W = 


E 

iei 


cPi 


a,- 


( 1 - 




■ ■ ■ (1 — X“ifc' 


This is a snfhciently general sitnation since we can always increase the nnmber of 
binomials in a fraction by nsing the identity 


xP _ xP(l-x“^) 

(1 — x“l) ■ ■ ■ (1 — x“'=-l) (1 — x“l) ■ ■ ■ (1 — x“'“) 

xP ^P+a-k 

(1 — X“l ) ■ ■ ■ (1 — X“'“ ) (1 — X“l ) ■ ■ ■ (1 — X“'= ) 


The procednre may increase the nnmber of terms by a factor of 2^, bnt since k is 
assnmed to be hxed, this amonnts to a constant factor increase. 

As nsnal, the inpnt size of (2.1) is the nnmber of bits needed to write /(x) down. 
Let /i,... , /d G E"" be integer vectors, U = (A^i,... , \in)- The vectors dehne the 
monomial map (f): C'^ as follows: 


( 2 , 2 ) 


(Cij-.-Xn)'—* ,a:d), whei-e x. ^ z' 


The inpnt size of this monomial map is the nnmber of bits needed to dehne it, that 
is, ronghly, dn + log 2 \Xij\. 

Snppose that the image of (f) does not consist entirely of poles of /(x). Then we 
can dehne a rational fnnction g : C by 

^(z) = /((/)(z)). 


The goal of this section is to constrnct a polynomial time algorithm, which, given 
a rational fu n ction (2.1) with a hxed nnmber k of binomials in each fraction and 
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a monomial substitution (2.2), computes a formula for g{z). Note that we cannot 
just substitute x = in the formula (2.1) since any z G may turn out to be 
a pole for some fraction of (2.1) and yet a regular point of g. For example, if d = 1, 
n = 0 and 


fix) 


X 


n+l 


1 — X 1 — X 


E 

Tn=0 


X 


then a: = 1 is the pole of both fractions but is a regular point of /; we have 

/(I) = n + 1. 

To this end, let us associate with the rational function (2.1) a meromorphic 
function F(c), c G C'^, dehned by 


(2.3) 


r(c) = 

iei 


_ exp{c,pi) _ 

(1 - exp(c, Uii)) ■ ■ • (1 - exp(c, atk))' 


As usual, for c G with c = r + it, where r, t G and a G we let (c, a) = 
(r, a)+i{t, a), where (■, ■) is the standard scalar product in The set of poles of the 
i-th fraction is the union over l<j<k of the hyperplanes {c G : (c, a^j) = O}. 
However, the set of poles of F{c) may be much smaller because of cancellations of 
singularities. 

There is a simple relation between (2.1) and (2.3). For c = (71 ,... , 7 ^) and 
X = (xi,... , Xd) we write 


X = provided Xi = exp{7i} for i = 1 ,... , d. 


Then the functions (2.1) and (2.3) are related by the equation 

rw = /(e“). 

Let L C be a subspace such that a generic c G L is a regular point of F{c). 
We want to construct a short formula for F{c) for c & L. We assume that the 
subspace L C is given by its integer basis. Again, we cannot just use (2.3), since 
L may be orthogonal to some vectors aij and hence a generic c & L may be a pole 
of some fractions in (2.3) while being a regular point of F{c). 


(2.4) Definition. Given /, let us consider the function 


.&) = nTV 




1 “ exp{-T^J 

in / + 1 (complex) variables r and ^1 ,... , It is easy to see that G is analytic in 
a neighborhood of the origin r = = ... = ^; = 0 and therefore there exists an 

expansion 


+ 00 


Gir-,^1,... ,6) = X]rWdj(^i,... ,^z). 


i=o 

where tdj(^i,... , ^z) is a homogeneous polynomial of degree j, called the j-th Todd 
polynomial in ^1 ,... It is easy to check that tdj(^i,... ,^z) is a symmetric 
polynomial with rational coefficients, cf. [BP99]. 
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(2.5) Lemma. Let us fix k. Then there exists a polynomial time algorithm, which, 
given a function (2.3) and a subspace L C which does not lie entirely in the set 
of poles of F computes F{c) for c E L in the form 


F(c) = ^ ft 


exp(c,Qi) 


iei' 


(l - exp(c, bii)) ■ ■ ■ (l - exp(c, bis)) ’ 


where s < k, fii E Q, qi, bij E and bij is not orthogonal to L for any i,j. 

Proof. Let us consider the representation (2.3). Let us choose a vector n G such 
that {v, aij) 7^ 0 for all Uij. Such a vector v can be constructed in polynomial time, 
see, for example, [BP99]. Let r be a complex parameter. Then, for any regular 
point c of F{c) the function F{c + rv) is an analytic function in a neighborhood of 
r = 0 and the constant term of its expansion at r = 0 is equal to F{c). Hence our 
goal is to compute the constant term (in r) of every fraction in the representation 
(2.3) of F{c + tv) and add them up. 

Let us consider a typical fraction 

^_ exp(c + TU,p) _ 

(l - exp(c + TV, ai)) ■ ■ ■ (l - exp(c + tv, a^)) ’ 

where p, aj E as a function of r. Suppose that the vectors ai orthogonal to L 
are ai,... ,ai for some I < k. Then 


h(r) 


i 

=T~^ exp(c, p) exp{T(u, p)} JJ 

i=l 


T 

1 - exp{r(u,ai)} 


k 

i=l+l 


1 

exp(c + TV, ai) 


Now we observe that fihij) is an analytic function of r and that our goal is to 
compute the coefficient of in the expansion of fihij) in the neighborhood of 
r = 0. 

First, we observe that 


+ 00 


(2.5.1) 


exp{T(u,p)} = ^ 


i=o 


{v,p)^ 


■tF 


Second, letting f^i = —{v,ai) for f = 1 ,... ,1, we observe that 


+ 00 


(2.5.2) 


n 1 _ exp{r(., a,>} ft • • ■ ft ''' ' 



Finally, 


(2.5.3) 


k 


n 

i=l+l 


1 

1 — exp(c + TV, tti) 


+ 00 

^ ^ Hj (c, , 

i=o 


,ak,v)T^ 


for some functions Hj. 

Note that r = 0 is a regular point of 


k 


n 

i=l+l 


1 

1 — exp(c + TV, tti) 


and so we compute Hj differentiating the product j times and setting r = 0. By 
the repeated application of the chain rule, Hj is a polynomial in exp(c, a^), {v, ai) 
and (1 — exp(c, ai))~^. Thus, for all ji, j 2 , Js such that ji + j 2 + js = I, we have to 
combine the ji-st term of (2.5.1) , the j 2 -nd term of (2.5.2) and the js-rd term of 
(2.5.3). Since I < k and k is hxed, we get the desired result. □ 

Remark. If L = {0} and 0 is a regular point of F[c), the algorithm of Lemma 2.5 
computes the number F{0). This procedure is used in [B94] to compute the number 
of integer points in a polytope. 

Now we can compute the result of the monomial substitution (2.2) into the 
rational function ( 2 . 1 ). 

(2.6) Theorem. Let us fix k. Then there exists a polynomial time algorithm, 
which, given a function (2.1) and a monomial map : C” —given by (2.2), 
such that the image of does not lie entirely in the set of poles o//(x) computes 
the function g{z) = /((/)(z)) as 


where s < k, fii E Q, qi, bij G IF and bij 7 ^ 0 for all i,j. 

Proof. Let F{c) be the function (2.3) associated to /(x). With the monomial map 
( 2 . 2 ) we associate a linear transformation $ : 

Cl—^ ((c,/i),... ,{c,ld)) 

and the adjoint transformation : 

^*(^1, • ■ ■ , ^d) = ^ 1/1 + • • • + Cdld- 

G{c) = F($(c)) for c e C^. 
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Let us dehne 



Hence 


G'(c) = g{e^). 


Let L C C'^ be the image of under $. Then L does not he entirely in the set of 
poles of F{c). Applying Lemma 2.5, we compute G{c) = F($(c)) in the form 


G(c) = 5^/3. 
iei' 


_ exp($(c),'u^) _ 

(1 - exp($(c), Vii)) ■■■{!- exp($(c), Vis)) ’ 


where for i,j we have {^{c),Vij) ^ 0 for a generic c & L. Now we let qi = 
and bij = ^*{vij) so that 




G(c) = 5^ ft 
iei' 


_ exp(c, Qi) _ 

(1 - exp(c, 6ii)) ■ ■ ■ (1 - exp(c, bis)) 


and the result follows. □ 

Remark. In particular, if x = (1,... , 1) is a regular point of (2.1), we can choose 
/i = ■ ■ ■ = /d = 0 in (2.2). In this case, the algorithm of Theorem 2.6 computes the 
value of /(I,... ,1). 

3. Operations with Generating Functions 

Some of the results of this section are stated in [BP99]. Many of the proofs in 
[BP99] are only sketched and some non-trivial details are omitted. We give a mostly 
independent presentation with complete proofs. The main goal of this section is to 
prove that if hnite sets 81,82 C are given by their generating functions f{ 8 i-, x) 
and /(S' 2 ;x) then the generating function f{ 8 ;yi) of their intersection 8 = 81 H 82 
can be computed efficiently. Our main tool is the generating function for the integer 
points in a rational polyhedron. 

Let P C be a rational polyhedron and let S' = P fl be the set of integer 
points in P. Let 

/(S;x)= ^ x™, 

mePnzd- 

Thus if P is bounded, /(S;x) is a Laurent polynomial in x. If P (possibly un¬ 
bounded) does not contain straight lines then there is a non-empty open set U G 
such that the series converges absolutely and uniformly on compact subsets of U 
to a rational function of x. If P contains a straight line it is convenient to agree 
that /(S;x) = 0, see [BP99]. 

We need the following result from [BP99] which states that /(A; x) can be written 
as a short rational function. 
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(3.1) Theorem. Let us fix d. Then there exists a polynomial time algorithm, 
which, for any given rational polyhedron P dMfi computes f{P fl x) as 


/(Pnz‘';x) = 


cPi 


iei 


(1 — x“»l ) ■ ■ ■ (1 — x“i<i) 


where Ci G {—1,1}, Pi,aij,& Tfi and aij 7 ^ 0 for all i,j. In fact, for each i, 
an,... , aid is a basis ofh^. 

A (complete) proof can be found in [BP99], Theorem 4.4. 

To compute the generating function of the intersection of two sets, we compute 
a more general operation, that is, the Hadamard product of two rational generating 
functions. 

(3.2) Definition. Let gi and g 2 be Laurent power series in x G 


5 fl(x) = ^ AmX"" and 5 f 2 (x) = ^ /32mX"'. 


The Hadamard product g = gi* g 2 is the power series 

5 f(x) = ^ /3mX"^ where (3rn = /3lm(32m 

maT/- 


First we will show that the Hadamard product of the Laurent expansions of 
some particular rational functions can be computed in polynomial time. Namely, 


let us choose a non-zero vector I G 


and suppose that an,... ,aik G Z“ and 


« 2 i, • ■ ■ , ci 2 k £ are vectors such that (/, a^j) < 0 for all i,j. Let pi,P 2 G and 
let 


(3.3) ^i(x) = 


rPl 


(1 — X“ll) ■ • • (1 — x^ifc) 


and g 2 (x) = 


rP2 


(1 — X“2l) ... (1 — x“2fc) 


We observe that for all x in a sufficiently small neighborhood t/ of xq = e^, we have 
|x“o | < 1 and so gi and g 2 have Laurent series expansions for x G 17. Indeed, if 
|x“| < 1 , the fraction 1/(1 — x“) expands as a geometric series 


1 


1 — x“ 




and to obtain the expansions of g\ and g 2 we multiply the corresponding series. 
Clearly, the Hadamard product of the expansions converges for all x G 17 to some 
analytic f un ction h, which we also denote gi* g 2 - We prove that once the number 
k of binomials in (3.3) is fixed, there is a polynomial time algorithm for computing 
the Laurent expansion of /i = ★ 5 f 2 as a short rational function. 
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(3.4) Lemma. Let us fix k. Then there exists a polynomial time algorithm, which, 
given functions (3.3) such that for some I & we have {aijfi) < 0 for all i,j, 
computes a function /i(x) in the form 


K^) = 

iei 


(1 — X^il) ■ ■ ■ (1 — X^i») 


with qi,bij G G Q and s < 2k such that h has the Laurent expansion in a 

neighborhood U o/xq = and /i(x) = ( 7 i(x) -k g 2 (x). 

Proof. In the space = |(^i,... ,^ 2 fc)} let P be a rational polyhedron dehned 
by the eqnations 


Pi + + • • • + ikO'lk — P2+ ik+ia'21 + • • • + i2k(^2k 

and the ineqnalities 

> 0 for i — 1,... , 2k. 

Let z = {zi,... , Z 2 k) and let ns consider the series 

(3.4.1) /(PnZ2'=;z)= 

mePnz^fe 

Clearly, the series converges absolntely and nniformly on compact sets as long as 
\zi\ < 1 for i — 1,... , 2k. By Theorem 3.1 we compnte f{P fl z) in the form 

(3.4.2) f(P n Z^'^; z) = V Ci--^^^ 

^ ^ ^ ^ ... (1 _ 

for some vectors Ui, Vij G Z^^ and some nnmbers G {—1,1}, where Vij 7 ^ 0 for all 
hi- 

On the other hand, expanding (71 (x) and (72 (x) as prodncts of geometric series, 
we obtain 


k 

^l(x) =xPl JJ ^ + 

j=lrti€Z+ (^1,... 

k 

^ 2 (x) =xP 2 JJ ^ = Y 

i=li^i€Z+ {v-i,... ,vu)el2^ 


Since the Hadamard prodnct is bilinear and since 


x^i 


★ X 


m 2 


X™.1 if 


0 if mi 7 ^ m 2 , 
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we conclude that 


(m,n)ePnZ^'' 

n={i/i,... ,Uk) 

Thus /i(x) is obtained from the function f{P fl Z^^;z) (cf. (3.4.1)-(3.4.2)) by 
the monomial substitution 

= x“ii ,Zk= x“i'=, Zk+1 = 1, . . . , Z2A: = 1- 

Now we use Theorem 2.6 to compute the result of the monomial substitution in 
(3.4.2). □ 

Now we are ready to prove the main result of this section. Suppose we have two 
hnite sets Si, S 2 C and let f{Si; x) and /(S'2; x) be the corresponding generating 
functions 

/(S,;x)= x™ and /(Sj^x) = x™. 

tuESi mES2 

Suppose further, that /(S'i;x) and /(S'2;x) can be written as short rational func¬ 
tions 


(3.5) 


/(Sdx) = 


i-pi 


ieh 


(1 — X“il) ■ ■ ■ (1 — X' 


^ik 


/(S,;x) = ^ft—^ 




iEl2 


(1 — X^»l) ■ ■ ■ (1 — X^ 


and 


with cti, (3i G Q, Pi, Qi, ttij, bij G Z'^ and aij, bij 7^ 0. Now let us consider Si and S 2 
as defined hj representations (3.5) of /(S'i;x) and /(S'2;x) as rational functions. 
Let S' = S'! n *52. Our goal is to compute the representation of 

/(S;x)= J^x™ 

mES 


as a short rational function. Again, we assume the number of k of binomials in each 
fraction of (3.5) hxed and allow numbers ai and fii and vectors pi, Qi and aij, bij to 
vary. 

( 3 . 6 ) Theorem. Let us fix k. Then there exists a polynomial time algorithm, 
which, given /i(Si;x) and /2(S2;x) computes f{S;x) for S = Si H S 2 in the form 


fiS;x) = 

iEl 


(1 


rVil 


)■■■{! 


0 ’ 
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where s < 2k, 7 ^ G Q, Ui, vij G and Vij 7 ^ 0 for all i,j. 

Proof. Let us choose a vector I G Z'^, such that (/, Oij) 7 ^ 0 and (/, bij) 7 ^ 0 for all 
i,j. As we remarked before, such a vector I can be constructed in polynomial time. 
When (/, ttij) > 0 or when (/, 6 ^ 7 ) > 0 we apply the identity 

xP _ xP-“ 

1 — x“ 1 — x““ ’ 

to reverse the direction of aij or bij, so that we achieve (/, aij) < 0 and (/, bij) < 0 
for all i,j in the representations (3.5). Then we can write 

f{Si; x) = ^ aifi'H(x) and f{S 2 ; x) = ^ Afi' 2 i(x) 
ieh ieh 

for some functions gii,g 2 i of type (3.3). There are Laurent series expansions of 
/(5'i;x) and /(5'2;x) in a neighborhood U of the point xq = and 

/(*5;x) =/(*Si;x) ★/(A 2 ;x) = ai^Pi^gii{x)-kg 2 iix). 

il€li,i2€l2 


We use Lemma 3.4 to compute f{S; x). □ 

Let Si,... , Sm C Z*^ be sets. We say that A C Z'^ is a Boolean combination of 
Ai,... , Sm provided S is obtained from Si by taking intersections, unions and com¬ 
plements. An immediate corollary of Theorem 3.6 is that the generating function 
of a Boolean combination of sets can be computed in polynomial time. 

(3.7) Corollary. Let us fix m (the number of sets Si C Z'^J and k (the num¬ 
ber of binomials in each fraction of f {Si] iL)). Then there exists an s = s{k,m) 
and a polynomial time algorithm, which, for any m (finite) sets Si,... , Sm C Z'^ 
given by their generating functions /(5'i;x) and a set S <Zl/' defined as a Boolean 
combination of Si ,... , Sm, computes /(5';x) in the form 


/(‘5;x) = 2^7^773 

iei 


(1 — x^*i) ■ ■ ■ (1 — x'^*®) ’ 


where 7 ^ G Q, Ui, Vij G Z'^ and Vij 7 ^ 0 for all i,j. 

Proof. We note that 

f{Si U S 2 ] x) = f{Si] x) f{S 2 ] x) - f{Si n S 2 ] x) and 
f{Si\S2]^) = f{Si]^)-f{SinS2]^) 

for any two subsets Si, S 2 C Z^. The proof follows by Theorem 3.6. □ 

Finally, we discuss how to patch together several generating functions into a 
single generating function. 
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(3.8) Definitions. By the interior int P of a polyhedron P C we always mean 
the relative interior, that is, the interior of P with respect to its affine hnll. 

Let X C be a set. We denote by [W] the indicator fnnction [W] : R, 




1 if rr e W 
0 x^X. 


We will need a simple formnla for the indicator of the relative interior of a polytope: 


(3.8.1) [intP] = (-1)‘^™^^(-1)‘^'”-^[P], 

F 


where the snm is taken over all faces of P inclnding P itself. This is a simple 
corollary of the Enler-Poincare formnla; see, for example. Section VI.3 of [B02]. 

From Theorem 3.1 we dednce the following corollary. 


(3.9) Corollary. Let us fix d. Then there exists a polynomial time algorithm, 
which, for any given rational polytope P C R*^ computes f{S; x) with S = (intP) fl 
Z'^ in the form 


/(‘5;x) = 


E 

iei 


rVi 


aj 


(1 — X“il) ■ ■ ■ (1 — X“i<i) ’ 


where ai G Q, pi, atj G if and aij 0 for all i,j. 
Proof. Applying formnla (3.8.1), we get 


/(S; x) = F/(F n Z-i; x). 

F 


Since the dimension d is fixed, there are polynomially many faces P and their 
descriptions can be compnted in polynomial time from the description of P. We 
nse Theorem 3.1 to complete the proof. □ 

Let ns consider the following sitnation. Let S' C Z*^ be a finite set and let 
Qi ,... , Qn C R^ be a collection of rational polytopes snch that S C ljr=i ™^Qi 
and int Qi fl int Qj = 0 for f 7^ j. In a typical sitnation, Qi,... , is a polytopal 
complex, that is, the intersection of every two polytopes Qi and Qj, if non-empty, 
is a common face of Qi and Qj and a face of a polytope Qi from the collection is 
also a polytope from the collection (in particnlar, not all Qi are fnll-dimensional). 
In this case, ur=l Qr = ljr=i Qi int Qi are pairwise disjoint. 

Snppose that we are given the fnnctions 


f{SnQj-,iL) 


E 

ieln 


i-Pi,. 


^i, 


(1 


rdil 


) . . . (1 — j 


and that we want to compnte /(S;x). In other words, we want to patch together 
several generating fnnctions /(S n(5j;x) into a single generating fnnction /(S;x). 
We obtain the following resnlt. 
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(3.10) Lemma. Let us fix k and d. Then there exists a polynomial time algorithm, 
which, given polytopes Qi,... , Qn and functions f{S fl <5j; x) computes /(S'; x) in 
the form 

/(S;x) = Vft- -- j— 

for s < 2k. 


Proof. We can write 


On the other hand, 


/(*5;x) = E f{S n int Qi] x). 

i=l 


S n int Qi = {S n Qi) n (int Qi fl Z^). 

First, nsing Corollary 3.9 we compnte f {int Qi fl Z^;x), and then nsing Theorem 
3.6 we compnte f{S fl int Qp, x). □ 


4. Lattice Width and Small Gaps 

In this section, we establish a simple geometric fact which plays a crncial role in 
the proof of Theorem 1.7. We start with dehnitions. 

(4.1) Definitions. Let A C be a lattice (that is, a discrete additive snbgronp 
of of rank d) and let A* C be the dnal (reciprocal) lattice, that is, 

A* = |c G M''* : {c,x) G Z for all x G a|, 

where (■, ■) is the standard scalar prodnct in Mfi. For a convex body B C (by 
which we mean a convex compact set) and a non-zero vector c G A* let 

width(i?, c) = max(c, x) — min(c, x) 

xEB xEB 

be the width of B in the direction of c. Let 

widthlS) = min widthlS, c) 
c6A*\{0} 

be the lattice width of B. As is well known, the minimnm indeed exists. 

It is known that there exists a constant Lv{d) with the following property: if 
B n A = ij) then width(i?) < u{d). It is conjectnred (and proved in many special 
cases) that u{d) = 0{d) while the best known valne is u}{d) = 0{d\nd) [BL-l-99]. 

We state some obvions properties of the width: 

width(i?, c) = width(i? -|- x, c) for any a: G M''* and 
width(ai?, c) = a width(i?, c) for all a > 0. 

Conseqnently, 

width(i?) = width(i? -|- x) for any a: G and 

width(ai?) = Q;width(i?) for all a > 0. 
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(4.2) Lemma. Let B C be a convex body, let c E be a non-zero vector and 
let 

7 min = min(c, x) and 7max = max(c, x). 

xEB xEB 

Let 7min < 7i < 72 < 7max bc numbcrs. Then there exists a point xq E B and a 
number 0 < a < 1 such that for 

A = a{B — xo) + xo = aB + (1 — a)xo 


one has A C B and 


min(c, x) = 7i and max(c, x) = 72. 

xEA xEA 


Proof. Translating B, if necessary, we can assnme that 7min = 0 . Dilating B, 
if necessary, we can assnme that 7max = 1 - Then 0 < 71/(1 — 72+71) < 1 , 
and, therefore, we can choose xq E B snch that (c, xq) = 7i/(l — 72 + 7 i)- Let 
a = (72 — 71). Then, for A = a{B — xq) xq = aB + (1 — a)xo, we have 


min(c, x) 

xEA 


(1 - a)7i 

1-72 + 71 


and 

/ \ , (l-a)7i 

max(c, X) — a -\ -= 72. 

xea 1-72 + 71 

Since B is convex, we have A G B. □ 

Now we can prove the main resnlt of this section. 

(4.3) Theorem. Let S C 6 e a convex body and let A C be a lattice. Let 
c E A* be a non-zero vector. Let us consider the map: 

(f) : B n A —Z, (f){x) — (c, x) 

and let Y = (f>{B fl A). Hence Y C Z is a finite set. 

Suppose that 

width(S,c) < 2width(S). 

Then for any yi,y 2 ^ Y such that y 2 — yi > 2u{d) there exists a y E Y such that 

yi <y < y 2 - 

Proof. Snppose that snch a point y does not exist. Let ns choose any 0 < e < 1/2 
and let 71 = yi + e and 72 = 1/2 — e- By Lemma 4.2 there exists xq E B and a 
nnmber a > 0 snch that for A = a{B — xq) + xq, A C B, we have 

min(c, x) — 71 and max(c, x) =72. 

xEA xEA 
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Then there is no integer in the interval [71,72] which is a valne of (c, x) for some 
a: G -B n A. Hence A fl A = 0 . Therefore, we mnst have 

width(A) < u{d). 

On the other hand, since A is a homothetic image of B we have 

width(A) = a width(B) and width(A, c) = a width(B, c). 


Therefore, 

72 — 7 i = width(A, c) < 2 width(A) < 2u{d). 

Hence 1/2 — yi — 2e < 2u{d) for any e > 0 and y 2 — Vi < 2u{d), which is a 
contradiction. □ 

In other words, the set Y G Z does not have “gaps” larger than 2u){d). We will 
nse the following corollary of Theorem 4.3. 

(4.4) Corollary. Let Y C Z be the set of Theorem f.3 and let m = \2u{d)]. For 
a positive integer I, let Y + I = \^y + I : y G T} denote the translation ofY by 1. 
IfY^th then the set 

m 

Z = Y\{j(Y + l) 

1 = 1 

consists of a single point. 

Proof. By Theorem 4.3, we have Z = {z}, where z = min{i/ : y gY}. □ 

5. Projections and Partitions 

In this section, we snpply the remaining ingredient of the proof of Theorem 1.7. 
This ingredient, np to a change of the coordinates, is a weak form of a lemma of 
R. Kannan [K92]. 

We describe it below. Let T : be a linear transformation snch that 

T(M‘^) = and T(Z'^) C Z^. Thns k < d and the matrix of T is integral with 
respect to the standard bases of and Then kerT is a rational (d — k)- 
dimensional snbspace of (that is, a snbspace spanned by integer vectors) and 
A = Z*^ n (kerT) is a lattice in kerT. As is known (see, for example. Chapter 1 of 
[C97]), a basis of A can be extended to a basis of Z*^ and hence any linear fnnctional 
i : kerT —M snch that ^(A) C Z can be represented in the form i{x) = (c, x) for 
some c G Z'^. The representation, of conrse, is not nniqne as long as kerT 7 ^ M'^. 
For c G (kerT)-*- (the orthogonal complement of kerT), the corresponding linear 
fnnctional is identically 0 . 

Let T C be a rational polytope. For y G let ns consider the hber 

Ty = jx G T : T{x) = yj 
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of X. For c G Z'^ \ (kerT)-*- we define the width of Py in the direction of c as 

width(Py, c) = niax(c, x) — niin(c, x) 

XEPy XePy 

and we dehne the lattice width of Py as 

width(P„) = min widthiPj,, c). 

^ ceZd\(kerT)J- ^ 

We observe that the lattice width of Py so dehned coincides with the width (as 
dehned in Section 4), with respect to A of a translation Py C kerT. 

We need the following resnlt, which is a (rephrased) weaker version of Lemma 
3.1 from [K92]. It asserts, ronghly, that one can dissect the image T{P) into 
polynomially many (in the inpnt size of P and T) polyhedral pieces Qi and hnd for 
every piece Qi a lattice direction Wi snch that for all y & Qi the lattice width of Py 
is almost attained at Wi. 

(5.1) Lemma. Let us fix d. Then there exists a polynomial time algorithm, which, 
for any rational polytope T C and any linear transformation T : R^ such 

that T(R'^) = R^ and T{7fi) = , constructs rational polytopes Qi,... ,Qn C Mfi 

and vectors wi,... , Wn G Z^ \ (kerT)-*- such that 

(1) For each i — 1,... ,n and every y G Qi, 

either width(Ty, rc^) <1 or width(Ty, w^) < 2width(Ty); 

(2) The interiors int Qi are pairwise disjoint and 

n 

IJintg, =T(T). 

i=l 

Proof. Let ns constrnct a rational snbspace T C R'^ snch that V fl (kerT) = {0} 
and (kerT) + V = R'^. Then the restriction of T onto V is invertible and we can 
compnte a matrix L of the linear transformation R^ —V which is the right inverse 
of T. 

Snppose that the polytope P is dehned by a system of linear ineqnalities 

T = jx G R^ : Ax < fej, 

where A is an n x d integer matrix and b is an integer n-vector. Then the translation 
Py C ker T of Py is dehned by the system of linear ineqnalities 

Py = ^x e kerT : Ax < b — ATyj. 

As y ranges over Q = T{P), vector b' = b — ALy ranges over the rational polytope 
Q' = b — AL[Q) with dimQ' < k. Since width(Ty, c) = width(Ty, c) for all y G Q 
and all c and width(Ty) = width(T^), the resnlt follows by Part 3 of Lemma 3.1 of 
[K92]. □ 
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6. Proofs 


Now we are ready to prove Theorem 1.7. 

Proof of Theorem 1.7. Without loss of generality, we assume that T(M^) = 
Indeed, if im(T) ^ we consider the restriction T : im(T). After a 

change of the coordinates, the lattice A = Z^nim(T) is identihed with the standard 
integer lattice. 

The proof is by induction on dim(kerT) = d — k. 

Suppose that /c = d, so dim(kerT) = 0 and T : l/ —is an embedding. 
Let Cl,... , Cd be the standard basis of Z*^ and let ti = T(ei). Then /(S';x) is 
obtained from f(P fl Z'^;y) by the monomial substitution yi = x** and we use 
Theorems 2.6 and 3.1 to complete the proof. 

Suppose that d > k, so dim(kerT) > 0. Let Qi, ... , Qn C R^ be the polytopes 
constructed in Lemma 5.1. It suffices to compute the functions f{S fl (5i;x) for 
i — 1,... ,n and then, using Lemma 3.10, we can patch them together and obtain 
/('5;x). 

Let us consider a particular polytope Q = Qi and the corresponding intersection 
S' n Q. Let w = Wi, w E \ (kerT)-*- be a vector whose existence is claimed by 
Lemma 5.1. Let us consider the linear transformation 

f : R'^ —^ R^+i = R^ © R, f{x) = {T{x), {w, x)) 
and the projection 

pr : ^ Kt pr((i ,... ,ik+l) = (fl.... ,6). 

Finally, let P' = {x E P : T{x) E Q] and S = f{P' n Z'^) C 

Clearly, SflQ = pr{S) and dim(kerT) = d —/c —1, so we can apply the induction 
hypothesis to T and compute /(S;z), where z = (x,a:fc_|_i), Xk+i E C. Our goal is 
to compute /(Sfl Q; x) from /(S; z). To do that, we construct a subset Z C S such 
that the projection pr : Z —S fl Q is one-to-one, and then we obtain /(S fl Q; x) 
from f{Z; z) by substituting Xk+i = 1. 

For a positive integer /, let S + / denote the translation of S by / along the last 
coordinate. 


S + l = {(6, ....&. 6+1 + 1): Ki,. .. ,6) € s}. 

Clearly, 

/(S + i;z) = l'n.i/(S;z), 

Let m = [2c<;(d — k)~\ (see Section 4) and let us dehne 

m 

Z = S\[J{S + l). 

Z=1 
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Using Corollary 3.7, we compute f{Z;z). 

Now we claim that the projection pr : Z —S' fl Q is one-to-one. Let us consider 
the projection pr : S —S' fl Q. For a y G S let us consider the preimage Sy C S 
of y. We observe that 


*^ 2 /= {( 2/5 {w,x)) : xePyDZ'^Y 

that is, Sy consists of all pairs (y, (rc,^)), where x is an integer point from the 
hber Py of P over y: 

Py = jx e P : T{x)=yy 

By Lemma 5.1, we have either width(Py, rc) < 1 or width(Py, rc) < 2width(Py). If 
width(Py, w) < 2width(Py), then, by Corollary 4.4, the set 

m 

Zy=S,\[j(S,+l) 

1=1 


consists of a single point, that is, the point of Sy with the smallest last coordinate. 
If width(Py, rc) < 1 then Sy consists of a single point and so Zy consists of a single 
point as well. Thus, in any case, for any y ^ Sr\Q the preimage Zy of the projection 
pr : Z —> S' n Q consists of a single point, so pr : Z —S fl Q is indeed one-to- 
one. Hence, using Theorem 2.6, we compute /(S fl Q;x) by specializing f{Z;z) at 
Xfc+i = 1 (where z = (x,a:fc+i)). □ 

We deduce Theorem 1.5 from Theorem 1.7. 

Proof of Theorem 1.5. Let us dehne a linear transformation T : R by 

... , = ai^i -F ... -F ttd^d- 

Thus S = is the semigroup generated by ai,... , a^. It remains to notice 

that there are some explicit bounds for the largest positive integer not in S, so one 
can replace the non-negative orthant by a rational polytope to get the initial 
interval of S. For example, in [EG72] it is shown that if t > max{ai,... , ad} then 
all numbers greater than or equal to 2t‘^/d are in S. Let n = \2t‘^/<f\ and let 

d 

P = |(^i,... <n-l and > 0 for f = l,...d| 

i=l 

be the simplex in R^. Then we can represent S as a disjoint union of T[P fl Z'^) 
and the integer points in the ray [n, -|-cx)) Since the generating function of the set 
of integer points in the ray [n, -|-cx)) is just x ^^^/(1 — x), applying Theorem 1.7 we 
complete the proof. □ 
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7. Further Examples: Hilbert Bases, Test Sets and Hilbert Series 

As another application of Theorem 1.7, let ns show that certain Hilbert bases are 
ennmerated by short rational fnnctions. 

Let Ri,... , Rd C be linearly independent vectors and let 


a 

H = ctiRi : 0 < cti < 1 for i = 1,... , dj 


i=l 


be the parallelepiped spanned by ri, ... , r^, and let K be the convex cone spanned 
by Ri,... ,Rd: 


d 

K = : cti > 0 for i = 

i=l 

We say that a point n G H fl is indecomposable provided v cannot be written in 
the form v = Vi + V 2 , where Vi and V 2 are non-zero integer points from H. The 
set S of all non-decomposable integer vectors in H is called the (minimal) Hilbert 
basis of the semigronp AT fl Z'^, since every integer vector in K can be written as a 
non-negative integer combination of points from S', see Section 16.4 of [Sc86]. Let 
ns show that as long as the dimension d is hxed, the set S has a short rational 
generating fnnction. 

(7.1) Theorem. Let us fix d. Then there exists a number s — s{d) and a poly¬ 
nomial time algorithm, which, given linearly independent vectors ri ,... ,Ud G Z'^ 
computes the generating function /(S;x) of the (minimal) Hilbert basis S of the 
semigroup of integer points in the cone spanned by ui,... ,Ud in the form: 

ai- -7-^-;- r—, 

{1 -■ ■ ■ {I - 

where I is a set of indices, ai are rational numbers, pi, bij G Z'^ and bij y 0 for all 
hi- 

Proof. Let ns constrnct a rational polyhedron Q G H which contains all integer 
points in H except 0. This can be done, for example, as follows: we constrnct vectors 
/i, ... , /d G Z*^ snch that (Z^, Uj) = 0 for i y j and (/j, r^) > 0, let / = /i Id 

and intersect H with the halfspace {l,x) > 1. 

Let P = QxQcM.‘^®Mfi = and let T : P —> be the transformation, 
T{x,y) — X -\- y. Let Si = T{P fl 7?^ and let S 2 = Q fl Z*^. Then the minimal 
Hilbert base S can be written as S = S 2 \ Si. The proof now follows from Theorem 
1.7 and Corollary 3.7. □ 

Yet another interesting class of sets having short rational generating fnnctions 
is that of “test sets” with respect to a given integer matrix. 
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(7.2) Test Sets. Let us choose anx d integer matrix A such that for any b G 
the polyhedron 

Pb = jx G : Ax < fej, 

is bounded. A point a G is called a neighbor of 0 with respect to A provided 
there is a polytope Pj, containing 0 and a and not containing any other integer 
point in its interior. The set 5'(A) of all neighbors of the origin is often called a 
test set. Test sets 5'(A) play an important role in parametric integer programming 
[S97]. The set S{A) is hnite, and it has some interesting (for d > 2) and not quite 
understood (for d> 3) structure. One can show that for any hxed d and n, given A, 
the generating function /(A; x) for S = S{A) can be computed in polynomial time 
as a short rational function. The proof follows from Theorem 7.1 and Corollary 3.7 
in a similar way as above, since S can be expressed as a Boolean combination of 
projections of sets of integer points in some rational polytopes. 

We note that other types of test sets studied in the literature, such as Schrijver’s 
universal test set and Graver’s test set (see [T95] and [St96]), also admit a short 
rational generating function. 

Finally, we describe one related problem of computational commutative algebra. 

(7.3) Hilbert series of rings generated by monomials. Let us consider integer 
vectors ai,... , G with non-negative coordinates and let S be the semigroup 
generated by ai,... , a^: 


d 

S ^ ^ • fJ'i G Z-|-^ . 

i=l 

Thus S can be represented as the image T(Z^^) under the linear transformation 

T : ^ M", T(ei,... , + • • • + 

The generating function /(A; x) can be interpreted as the Hilbert series of the Z^- 
graded ring R = C[x“L ... ,x“‘^], cf. [BS98] and Chapter 10 of [St96]. The set S is 
inhnite and Theorem 1.7 is not directly applicable (although it allows us to claim 
the intersection of S with any given polytopal region has a short rational generating 
function). However, one can still compute the whole function /(S'; x) in polynomial 
time as a short rational function provided the dimension k and the number d of 
generators are hxed. We also note that by applying a monomial specialization of 
f{S',x.) we can obtain the Hilbert series of R under a coarser grading. 

We sketch an algorithm for computing /(S;x) below. 

Without loss of generality we assume that 7 ^ 0 for i = 1,... , d. Let us consider 
the product 


^(S;x) = /(S;x)(l-x“^)---(l-x“‘^). 
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It is not hard to prove that 5r(5';x) is, in fact, a polynomial in x. This follows, for 
example, from the interpretation of /(5';x) as a Hilbert series, cf. Section 1.9 of 
[E95]. 

We need to compnte a bonnd L with the property that that if the coefficient of 
x"^, m = (/Ui,... in fi'(‘S';x) is non-zero then ni + ... + Hk < L. Snppose for 

a moment that we can hnd snch an L. Let ns consider the integer cnbe 

CCZ^, (7= {(^ 1 ,... ,ed) : 0<Cr<L for f = 1,... , d} 

and the integer simplex 

A C : A = |(//i, .. . ,//fc) : Hi + . . . + Hk < L'^■ 

Let S' = T(C'), so A n S' C S". Applying Theorem 1.7, we compnte /(S';x) as a 
short rational fnnction. Let 

< 7 (S';x)=/(S';x)(l-x“i)---(l-x“^). 


We note that 

d(‘5;x) = ^(S';x) ★/(A;x) 

and nse Theorem 3.1 and Lemma 3.4 to compnte the Hadamard prodnct g{S;yi) 
as a short rational fnnction. Finally, we let 


/(S;x) = ( 7 (S;x) 

7=1 

It remains, therefore, to compnte the bonnd L on the total degree of a monomial 
x"^ which may appear with a non-zero coefhcient in the expansion of g{S', x). 

Let ns consider the rational cone AT C 

iL=|(a:,y): and T{x)=T{y)^. 

The lattice semigronp K fl is hnitely generated and nsing some standard tech- 
niqnes (see Chapter 17 of [Sc86] and Chapter 4 of [St96]) one can compnte in poly¬ 
nomial time an npper bonnd M on the coordinates of generators yi) of ATflZ^'^. 
Let A be the snm of the coordinates of ai,... , aa- We claim that L = A{M + 1) is 
the desired npper bonnd. 

Indeed, for every generator (xi, yi) with Xi yi, let Zi — Xi — yi or Zi = yi — Xi, 
whichever is lexicographically positive. Thns each coordinate of Zi does not exceed 
M. Let 

z = zi\lj(zi + z,). 

i 
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One can observe that the restriction T : Z — S is one-to-one. In fact, for every 
X & S the vector z & Z snch that T{z) = x is the lexicographic minimnm among 
all y e snch that T{y) = x. 

For / C {1,... ,d} let C be the coordinate semigronp consisting of the 
points (^ 1 ,... , Cd) snch that = 0 for f ^ I. As is proved in [Kh95], the set Z can 
be represented as a hnite disjoint nnion of sets Zj of the type Vj + so that the 
coordinates of Vj do not exceed M. Let Sj = T{Zj). Then S is the disjoint nnion 



The snm of the coordinates of T{vj) does not exceed MA. Therefore, if x”^, m — 
(//i,... , Hk) appears with a non-zero coefRcient in the prodnct 




we mnst have < Af A -|- A = L, which completes the proof. 
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